Circular geodesies in the Hartle-Thorne metric 



M. A. Abramowicz * ^[ § G.J.E. Almergren * f 
W. Kluzniak * J A.V. Thampan f 

If Department of Astronomy and Astrophysics, Chalmers University, 

S-412 96 Goteborg, Sweden 
§ Silesian University at Opava, Department of Physiscs, 

CZ- 74601 Opava, Czech Republic 
f International School for Advanced Studies (ISAS/SISSA), 

Via Beirut 2-4, 1-34013 Trieste, Italy 
X Institute of Astronomy, University of Zielona Cora, 

Lubuska 2, PL-65-265 Zielona Cora, Poland 
JCopernicus Astronomical Center, ul. Bartycka 18, 00-716 Warszawa, Poland 
* UK Astrophysical Fluids Facility (UKAFF), Leicester University, England 

Abstract. The Hartle-Thorne metric is an exact solution of vacuum Einstein 
field equations that describes the exterior of any slowly and rigidly rotating, 
stationary and axially symmetric body. The metric is given with accuracy up 
to the second order terms in the body's angular momentum, and first order in 
its quadrupole moment. We give, with the same accuracy, analytic formulae 
for circular geodesies in the Hartle-Thorne metrics. They describe angular 
velocity, angular momentum, energy, epicyclic frequencies, shear, vorticity and 
Fermi-Walker precession. These quantities are relevant to several astrophysical 
phenomena, in particular to the observed high frequency, kilohertz Quasi Periodic 
Oscillations (kHz QPOs) in the X-ray luminosity from black hole and neutron star 
sources. It is believed that kHz QPO data may be used to test the strong field 
regime of Einstein's general relativity, and the physics of super-dense matter of 
which neutron stars are made of. 
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1. The astrophysical motivation 

Thanks to precise timing of binary pulsars' orbital decay pQ, the weak field limit of 
Albert Einstein's general relativity is far better tested than any other physical theory. 
However, the most interesting and bizarre predictions of Einstein's theory deal not 
with the weak field, but with the extremely strong field regime and these have never 
been tested. Thus, the very question |2] "Was Einstein right?" remains unanswered. 
One may argue that in the foreseeable future there will be no way to test super-strong 
gravity. Obviously, this cannot be done in laboratories on Earth. Central regions of 
Galactic sources containing black holes and neutron star do have gravity sufficiently 
strong for such tests, but they are only several tens of kilometer across and typically 
observed from kiloparsecs away. Thus, they cannot be spatially resolved with current 
instruments. However, already existing technology [2] provides very precise resolution 
in time of the observed variability of X-ray radiation coming from the vicinity of black 
holes and neutron stars. 

Of special importance here are the kilohertz, Quasi Periodic Oscillations (kHz 
QPOs) of the X-ray flux from galactic low mass X-ray binaries in which the compact 
object is either a black hole or a neutron star 0|. The QPOs are often observed in 
twin pairs of kilohertz frequencies. The observed frequency ratios are rational, with 
the particular ratio 3:2 being most common. In several sources it was possible to 
detect the orbital binary motion, and deduce the mass of the compact object from 
Kepler's laws. In a few cases the deduced mass is above the Chandrasekhar limit, 
suggesting a black hole. 

Rational ratios of oscillation frequencies suggest a resonance. There is a strong 
reason to think of a particular resonance here, the non-linear parametric resonance 
between the vertical and radial epicyclic oscillations in accretion disks 0- The world 
lines of matter in accretion disks differ only slightly from those corresponding to 
circular geodesies. The 3:2 resonance occurs at a precise location inside the disk, only 
a few gravitational radii from the source, i.e., in the super-strong gravity regime. This 
gives a unique opportunity for observational strong-gravity tests. In particular, for 
three black hole twin kHz QPO sources with known masses, the resonance model has 
already been used for direct measurements of the black hole spin [7] . That was 
possible due to the fact that the black hole metric is known analytically (the Kerr 
solution). 

For neutron and strange stars, the exterior metric is not analytically known in 
general. It is known only in the special but important case of a perfect fluid, stationary 
and axially symmetric star with mass Mq, angular momentum Jo, and quadrupole 
moment Qq, that rotates rigidly and slowly. In this case, the Hartle & Thorne [5] 
solution gives the metric analytically, with accuracy up to the second order terms 
in stellar dimensionless angular momentum j — J /AT 2 , and first order in terms of 
dimcnsionless quadrupole, q = —Q/M 3 . 

We give here, with the same accuracy, analytic formulae for circular geodesies 
in the Hartle-Thorne metrics. They describe angular velocity, angular momentum, 
energy, epicyclic frequencies, shear, vorticity and Fermi- Walker precession. These 
quantities are of interest in modeling of several astrophysical phenomena, in particular 
of the kHz QPOs (e.g., [§]). The radius of the innermost stable marginally orbit and 
related quantities have previously been computed to linear order in j by Kluzniak & 
Wagoner ^U], through quadratic and higher orders in j and q by Shibata & Sasaki [TT]. 
Sibgatullin & Sunayev [T2\, and numerically by Berti & Stergioulas |T5| . 
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Although the derivation of the following formulae is simple, it is technically 
troublesome and one must worry about typographical errors. We derived all the 
formulae in a semi automatic way |14) using Mathematica. 

2. The Hartle-Thorne Metric 

We use the geometrical units, in which c = 1 = G. In these units, mass, length 
and time are all measured in centimeters. Adding an asterisk (*) as the indicator of 
the usual physical units, we write for the mass, radius and time in geometrical units, 
M = GM*/c 2 , r — r* , and t — ct*. In these units and with spherical coordinates 
(t,r,9,(j>), the Hartle-Thorne metric |15l 1161 15| reads, 

ds 2 = g tt dt 2 + g rr dr 2 + g ee dd 2 + g<j><j> d(f> 2 + g<t>t d(f>dt + g^dtd0, ( 1 ) 

where g t x = g^ t and the components of the metric tensor are given below. 

gtt = +(l-2M/r) [l + j 2 F 1 + qF 2 
g„ = - (1 - 2M/r) 1 [1 + f Gx - qF 2 

- 2 [l+j 2 H 1 +qH 2 ] (2) 

St* = -2(M 2 /r)j S in 2 6 

where 

Fi = [8Mr 4 (r-2M)]~ 1 ■ 

[u 2 (48M 6 - 8M 5 r - 24A/V - 30M 3 r 3 - 60A/V + 135Mr 5 - 45r 6 ) 
+(r - M) (16M 5 + 8M 4 r - 10A/V - 30Mr 4 + 15r 5 )] + Ai(r) (3) 
F 2 = [8M r (r - 2Af )] _1 (5 (3u 2 - 1) (r - M) (2M 2 + 6Mr - 3r 2 )) - A x (r) (4) 
d = [8M r 4 (r - 2M)] ~ l {{L - 72M 5 r) - 3u 2 (L - 56M 5 r)) - A x (r) (5) 
L = (80A/ 6 + 8A/V + 10A/ 3 r 3 + 20Af V - 45Mr 5 + 15r 6 ) (6) 

15r(r-2Af)(l-3u 2 ) / r \ 
Al= 16M^ () 

i?i = (8Mr 4 )- 1 (1 - 3u 2 ) (16Af 5 + 8Af 4 r - lOAfV + 15A/r 4 + 15r 5 ) + A 2 (r) (8) 
H 2 = (8A.fr)" 1 (5 (1 - 3u 2 ) (2Af 2 - 3A/r - 3r 2 )) - A 2 (r) (9) 

15(r 2 -2M 2 )(3^ 2 -l) / r \ 
= 16AP hl J (10) 

and where in addition u := cos9. 

We have checked by directly calculating the Ricci tensor Rik, that the above 
Hartle-Thorne metric is indeed a solution of the vacuum Einstein's field equations: 
with the relevant accuracy to quadratic terms in j and linear terms in q, Rik = 0. 

The Kerr metric in the Boyer-Lindquist coordinates could be obtained from the 
above Hartle-Thorne metric after putting a — Mj, q = j , and making a coordinate 
transformation implicitly given by, 

r BL =r- a 2 /(2r 3 ) ((r + 2M) (r - M) + u 2 (r - 2M)(r + 3Af )) (11) 
6 BL =6- a 2 /(2r 3 ) (r + 2Af) cos (9 sin (9 (12) 
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3. The detailed formulae 

3.1. The Horizon 

The radius of the horizon (if present) is obtained by setting g 2 t — g tt g^ = and 
solving for r. 



2M 



l-^ 2 (7-15« 2 ) + ^<z(l-3u 2 ) 



(13) 



3.2. The Ergosphere 

Similarly, the radius of ergosphere can be easily obtained by setting g tt = 0, 



r Q = 2M 



(14) 



3. 3. Dragging of Inertial Frames 

The angular velocity of an inertial observer near the rotating star as observed by 
someone at infinity (i.e. the frame-dragging) is given by: 

g t4> _ 2 J _ 2M 2 j 

rtt 3 V* 3 



(15) 



As such, it is independent of the quadrupolc moment and the second order terms as 
well as the angle of inclination. 



3.4. The Orbital Angular Velocity (SI = vf/v?) 

The angular velocity for corotating/counterrotating circular particle orbits is given 
by: 



SI = ± 



M 1 / 2 

r 3/2 



M3/2 
r 3/2 



j 2 F 1 (r) + qF 2 (r) 



where 



Fi(r) = [48M 7 -80M 6 r + 4M 5 r 2 -18M 4 r 3 +40M 3 r 4 + 10M 2 r 5 
+ 15 Mr 6 - 15r 7 ] (16M 2 (r - 2M) r 4 )" 1 + A(r) 
5 (6 M 4 - 8 M 3 r - 2 M 2 r 2 - 3 Mr 3 + 3 r 4 ) 



_ 15 (r 3 -2M 3 ) 
A(rj _ 32 M3 



16 M 2 (r-2M) r 
r 



In 



2M 



5.5. The Specific Angular Momentum (I = —u^/ut) 



±£0 



1TJ 



. M 3 / 2 (3r-4M) , . . _ . 

+ .] 2 F 1 (r) -qF 2 (r) 



r 3/2( r _2M) 



(16) 

(17) 
(18) 

(19) 
(20) 

(21) 
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where 



t := 



Fi(r) 
A(r) 



M l/2 r 3/2 

r- 2M 

16M 2 r 4 (r - 2M) 2 ] 1 (96M 8 - W2M 1 r - 8M 6 r 2 - 48M 5 r 3 

+42MV + 220M 3 r 5 - 260M 2 r 6 + 105Mr 7 - 15r 8 ) + A{r) 

[16M 2 r (r - 2M)] _1 5 (6M 4 - 22M 2 r 2 + 15Mr 3 - 3r 4 ) + A(r) 
15 



32M 3 



(2M 3 + 4M 2 r - 4Mr 2 + r 3 ) In 



2M 



5.5. T/ie Specific Energy (e = u t ) 
e = E [1 T j F i (r) + j 2 F 2 (r) + g F 3 (r)] 
where 
So 



2M 



Fi(r) 
B(r) 



rV2( r _ 3M)V2 
rV2 ( r -2M)(r - 3M) 

[l6Mr 4 (r - 2M) (r - 3M) 2 ] _1 (144M 8 - 144M 7 r - 28M 6 r 
-58M 5 r 3 - 176MV + 685M 3 r 5 - 610M 2 r 6 + 225Mr 7 - 30r 8 ) + B(r) 
5 (r — M) (6M 3 - 20M 2 r - 21Mr 2 + 6r 3 ) 



16Mr(r-2M) (r - 3M) 
15r (8M 2 -7Mr + 2r 2 ) 



-B(r) 



32M 2 (r - 3M ) 



In 



r-2M 



(22) 

(23) 
(24) 

(25) 

(26) 

(27) 

(28) 
(29) 

I 

(30) 
(31) 



3. 7. Radius of marginally stable, marginally bound and photon orbit. 

The condition e = u t = 1 gives the radius of the marginally bound orbit, r m b and the 
condition u* = gives the photon orbit, r v h- In addition by setting dljdr = 0we can 
solve for the radius of the marginally stable orbit, r ms . 



r mb = AM [1 =f \ j - j 2 (8047/256 - 45 In 2) + q (1005/32 - 45 In 2)] 



r p h 



3M 



. 2 V3 



/ 7036- 6075 In 3 



1296 



= 6M 



i 1 

TJ 3V 3 +J V 2592 



2 / 251647 , 3 

2 ' - 240 In - 



7020- 6075 In 3 
1296 

9325 
96 



240 In 



(32) 
(33) 

(34) 



3.8. The Epicyclic Frequencies 

We define an effective potential as: U(r,0,£) := g a - 21 + l 2 g^, which can be 
used to find the general formula for the epicyclic frequencies on the equatorial plane. 



(g tt + tt ±gt0 ) 2 (d 2 U 



2g 2 



dx 2 



x e (r, 9) 



(35) 
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Then we have: 

k 2 = M(r - 6Ai>- 4 [l ± j F x (r) - j 2 F 2 {r) -qF 3 (r)] 
K 2 e = Mr- 3 [1 T j Gx (r) + j 2 G 2 (r) + g G 3 (r)] 
where 

6 Af 3 / 2 (r + 2M) 
FlW ~ r 3 / 2 (r - 6M) 

F 2 (r) = [8M 2 r 4 (r - 2M) (r - 6M)] _1 [384Af 8 - 720M 7 r - 112AfV - 76M £ 
- 138A/V - 130A/V + 635A/V - 375Mr 7 + 60r 8 ] + A(r) 

5 (48Af 5 + 30A/ 4 r + 26M 3 r 2 - 127Af 2 r 3 + 75Mr 4 - 12r 5 ) 

F-i(r) = — — A(r) 

3K ' 8M 2 r (r - 2M) (r - 6M) V ' 

15r (r - 2M) (2Af 2 + 13A.fr - 4r 2 ) , 
A(r) = y ^ — — In 



Gi(r) 



16Af 3 (r-6Af) \r-2M 
6 AT 3 / 2 



r 3/2 

G 2 (r) = [8M 2 r 4 (r - 2M)] ^ [48A/ 7 - 224M 6 r + 28Af 5 r 2 

+ 6Af 4 r 3 - 170A/ 3 r 4 + 295AfV - 165Mr 6 + 30r 7 ] - B(r) 

, . 5 (6Af 4 + 34Af 3 r - 59Af 2 r 2 + 33Mr 3 - 6r 4 ) 

^ = ^r(r-2Af) ~ + ^ 

15 (2r - Af) (r - 2Af ) 2 , / r 
B ( r ) = TaTH ln 



16M 3 V r - 2Af 

3.9. Shear and Vorticity 

The general formulae for shear and vorticity are |17|: 

a 1 = - \ (i-ney 2 n 2 (v a n)(v a n) 

oj 2 = -\{l-VLl)- 2 K.- 2 {V a l){V a l) 

\ E J g^-gttg^ 

The results for the HT metric are: 

a 2 = S Q [1 ± j fi(r) + j 2 F 2 (r) + q F 3 (r)] 

oj 2 = V [1 TjGi(r)+j 2 G 2 (r) + 9 G 3 (r)] 

where 

9M (r-2Af) 2 



5 :- 
Vo := 



16r 3 (r-3Af) 2 
M (r-6Af) 2 
16T 3 (r - 3M ) 2 
and 

. . . 4Af 3 / 2 
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F 2 (r) = [8M 2 r 4 (r - 2M ) 2 (r - 3Af) 2 ] -1 • 

- (864M 10 - 1776M 9 r + 1048MV - 592Af 7 r 3 - 10M 6 r 4 + 542M 5 r 5 

- 693M 4 r 6 + 820MV - 530M 2 r 8 + 150Mr 9 - 15r 10 ) - A(r) (54) 
, x 5(36M 6 - 2M 5 r - 28Af 4 r 2 + 35M 3 r 3 - 43Af 2 r 4 + 21Afr 5 - 3r 6 ) 4 , w . 

^ - J 8MV(,-2M) 2 (,-3M) " + A(r)(55) 

= 15 r(r -4M)(M 2 + , 2 ) 
W 16A/ 3 (r-3Af) V r - 2M / 

12M 3 / 2 (r - 2M) 

ltrj_ r 1 /2( r -3M)(r-6M) ( ' 

G 2 (r) = - (7776A/ 10 - 12528Af 9 r + 3672A/V - 4080MV + 150Af 6 r 4 

- 12594M 5 r 5 + 30891M 4 r 6 - 25620A/ 3 r 7 + 9590M 2 r 8 - 1650Mr 9 

+ 105r 10 ) • [8AfV (r - 3M) 2 (r - 6A/) 2 ] ~* + B(r) (58) 

^ , , -5(54M 5 + 30A/ 4 r + 147M 3 r 2 - 271Af 2 r 3 + 141Afr 4 - 21r 5 ) . . , . 

G3W = 8MMr-3M)(r-6M) W ( } 

15 r (24Af 4 - 126Af 3 r + 135A1 2 r 2 - 54Alr 3 + 7r 4 ) / r \ 

S(r) = 16A1 3 (r - 3M) (r - 6M ) ln ^r-2M J (60) 

$.10. TTie Fermi- Walker Precession 

The precession for a gyroscope on a general circular orbit in a general stationary and 
axially symmetric spacetime was calculated in [T^]. In the equatorial plane for the 
Hartle-Thorne metric we use formula (4.8) of to obtain: 

n = fl [-1T J F 1 (r)+j 2 F 2 (r)+qF 3 (r)] (61) 
M i/2 r 2 

11(5 (r-3Af)V2 ( r -2M) 3 ^ 

Fl{r) = 2r 3 /2 (r-2M) (r-3M) (63) 
F 2 (r) = (576M 10 - 1704Af 9 r + 2904Al 8 r 2 - 1280A/V + 1604M 6 r 4 - 8060Al 5 r 5 
+10688A/V 6 - 6290APV + 1865A/V - 270Mr 9 + 15r 10 ) • 
[16M 2 r 4 (r - 2Af ) 2 (r - 3M ) 2 ] _1 - A(r) (64) 

5 (12A/ 4 + 36M 3 r - 94A/ 2 r 2 + 33A/r 3 - 3r 4 ) 

= 16MMr-3M) + AW (65) 

15r (r-4Af)(10Af 2 -8Afr + r 2 ) l / r \ 

AW = 32M3 (r - 3M ) ^^^J (66) 

5.11. Resonance 



Assuming an m : n resonance between epicyclic frequencies, we calculate the radius 
for which the resonance occurs given the condition that mK r = nUg and denoting 
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2 / 2 

z := m /n , 



r [l - j Fx (m, n) + j 2 F 2 (m, n) + q F 3 (m, n 



)] 



(67) 



where 



r °= 6M ( 1 "3 



(68) 



and 




(69) 



F 2 = (1 - 50 z - 11 z 2 + 385 z 3 + 10612 z 4 + 123286 z 5 + 496927 z 6 



+691843 z 7 + 251647 z 8 ) • 1296 z 4 (z - l) 3 (1 + 2 z) 



A(z) 



(70) 



-5 (1 + 5 z) (l + 74 z + 546 z 2 + 950 z 3 + 373 z 4 




(71) 



15 (l + 2z) (l + 12z + 63z 2 + 32z 3 ) / 3 



(72) 



4. Software and hardware 

These results were all derived and checked using Mathematica 4.0 ^H] running on a 
PHI 800 MHz/256MB under Linux Red Hat 7.2. In addition the packages TTC [H] 
and GRTensorM [2] were found useful. The Mathematica notebooks containing the 
calculations can be found on: 
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